function z = findststNoL(mu)

global alpha betta depr sigmaC sigmaL omega ALPHA popweight N e_real
global phi govexp
global lamda Delta1

Delta2=-Delta1*phi(1)/phi(2:N);
rstst = 1/betta-1+depr; %r from Euler at the steady state
kstst=e_real/(rstst/alpha)^(1/(1-alpha)); % express k from r=F_k
wstst = (1-alpha)*kstst^alpha*e_real^(-alpha); %w=F_e
lstst=e_real;

cstst = 1/(popweight(1)+sum(popweight(2:N).*lamda))*(kstst^alpha*e_real^(1-alpha) - depr*kstst - govexp); %c from resource constraint

z = mu*(popweight(1)+sum(popweight(2:N).*lamda))-cstst^(-sigmaC)-sum(ALPHA.*lamda.*(lamda*cstst).^(-sigmaC)); %FOC for consumption
